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Abstract 

We extend the Vasicek loan portfolio model to a setting where lia- 
bilities fluctuate randomly and asset values may be subject to systemic 
jump risk. We derive the probability distribution of the percentage loss 
of a uniform portfolio and analyze its properties. We find that the impact 
of liability risk is ambiguous and depends on the correlation between the 
continuous aggregate factor and the asset-liability ratio as well as on the 
default intensity. We also find that systemic jump risk has a significant 
impact on the upper percentiles of the loss distribution and, therefore, on 
both the VaR-measure as well as on the expected shortfall. 
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1 Introduction 



Modeling correlation of defaults plays naturally a central role in the liter- 
ature on the management of loan portfolios and valuation of credit deriva- 
tives. The reason for focusing on modeling the correlation of defaults is 
obvious from a credit risk management perspective: since default prob- 
abilities are relatively low, the (tail) dependence of different risks may 
play a dominant role in the joint probability of defaults. The role of this 
dependence is naturally pronounced during financial crises and, therefore, 
incorporating a factor (or factors) taking into account such relatively rare 
phenomena is clearly of interest. 

In light of the fact that financial crises occur regularly and share sev- 
eral similar characteristics both in their precedents as well as in their 
negative impact on the overall economy (for excellent recent empirical 
studies on financial crises see Reinhart and Rogoff (2008a, 2008b, 2009)), 
we extend the classical Vasicek loan portfolio model to a setting where 
the liabilities are subject to random fluctuations and asset values may 
face systematic asymmetric jump risk. We follow a factorized approach 
(for a well-written introduction of this credit risk modeling approach, see 



Schonbucher (2001 1) and model the liabilities as geometric Brownian mo- 



tions driven by two continuous factors; one aggregate level factor affecting 
all assets and liabilities and one statistically independent idiosyncratic fac- 
tor. In this way liabilities become conditionally independent of each other 
given the aggregate factor dynamics. On the other hand, we assume that 
the value of the assets constitute exponential jump diffusions potentially 
driven by three factors. As in the case of liabilities, the continuous part 
of the fluctuations are determined by an aggregate and an idiosyncratic 
factor. In addition to these factors, the assets may also be subjected to a 
discontinuous and spectrally negative risk component modeled as a com- 
pensated Poisson process with only positive jumps. This risk factor is 
assumed to be systemic in the sense that it affects all assets and, once 



realized, may significantly decrease their values (see, for example, Eber- 
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lein and Madan ( 2009 1 for a recent study where assets and liabilities are 



modeled as exponential Levy processes). This assumption allows the anal- 
ysis of rare but potentially dramatic collapses in the values of the assets 
backing up liabilities. Instead of considering valuation issues, we focus on 
the determination of the probability distribution of the percentage loss of 
a uniform portfolio and investigate how the interplay between different 
risk components affect this distribution. 

It is worth emphasizing that our model is related to the pioneering 



work by |Zhou|fl997| ) and |Zhou| ( |200lj ). |Zhou| ( |T997| ) (see also |Zhou| ( |200Tj ); 

for a recent application of Zhou's model within consumer credit setting, see 



de Andrade and Thomas (2007 1 ) presents a factorized model of credit risk 



where the value of the shares of the firm are assumed to evolve according 
to a geometric jump diffusion and the threshold value increases at a known 
constant exponential rate. Thus, liability risk is not considered as in our 



study. Moreover, in both Zhou ( 19971 and Zhou (2001 1 the jumps in the 
value of shares may have both signs implying that even though the firm 
is subject to sudden large drops, it may also face unexpected significant 
increases in its value. In this way, his model does not take into account the 
realization of rare but potentially significant negative outcomes eroding 



the value of assets. However, in contrast to our study, Zhou ( 19971 derives 
also the arbitrage free bond prices both in the case where default occurs 
whenever the value falls short the known threshold value at a given fixed 
maturity as well as in the technically demanding first passage time setting 
where default occurs whenever the value falls short the threshold value 
prior expiry. He also studies the term structure of credit spreads and 
demonstrates that the considered class of models can generate a variety 
of yield spread curves as well as marginal default rate curves. 

Our findings indicate that the impact of liability risk on the distri- 
bution of the percentage loss is ambiguous and depends, among others, 
on the correlation between the continuous aggregate factor and the asset- 
liability ratio as well as on the default intensity. If liabilities are subject 
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to purely idiosyncratic risk and are unaffected by the aggregate market 
factor, then the resulting probability distribution is similar but yet not 
identical with the limiting distribution in the standard Vasicek setting. 
The main reason for this is that in the present setting also the idiosyn- 
cratic risk factor affects the total volatility of the percentage growth rate 
of the asset-liability ratio. However, in contrast with the standard Vasicek 
model, our results indicate that there are circumstances under which the 
percentage loss converges to the known default probability of an individ- 
ual loan. Such a case arises when the continuous aggregate factor affects 
equally strong both assets and liabilities. Such a case might potentially 
appear in situations where both assets and liabilities depend on a well 
diversified portfolio (for example, in the case of unit linked products). 
The impact of the systemic risk component is more pronounced and our 
findings show that its presence may have a radical impact on the limiting 
distribution. First, the limiting probability distribution may have more 
than two modes; a phenomenon which does not arise in the standard con- 
tinuous setting. Second, the systemic risk term has a significant impact 
on the tail probabilities and tends to increase both the upper percentiles 
as well as the expected shortfall associated to these percentiles even when 
the realization intensity is low. 

The contents of this study are as follows. In section two we present 
the basic continuous model and state our main findings on the probability 
distribution of the percentage loss of a uniform portfolio. In section three 
we introduce the discontinuous systemic risk component and analyze its 
impact on the probability distribution of the percentage loss. Finally, 
section four concludes our study. 

2 The Impact of Liability Risk 

Our main objective is to investigate the percentage loss distribution of 
a loan portfolio within a conditionally independent factor model along 
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the lines indicated by the pioneering work in Vasicek ( 1987 1 (see also 



Vasicek ( 1991 I and Vasicek (2002)). To this end, we first assume that the 



asset values evolve according to the random dynamics characterized by 
the stochastic differential equation 

dA it = fj,iA it dt + (TiAudWit, Am = At, i = 1, . . . , n, (1) 

where both the drift coefficient \ii as well as the volatility coefficient a 
are exogenously given and Wu is standard Brownian motion. In order to 
model the statistical dependence of the various asset values, we assume 
that the driving Brownian motions can be decomposed into the form 

Wu =JpiYt + y/l-pi X iu (2) 

where Yj, Xu, . . . , X nt are a family of independent driving Brownian mo- 
tions and pi G [0,1], i = l,...,n measures the correlation between the 
underlying driving factor dynamics. The factor Y t is a, joint aggregate risk 
factor (market risk) affecting all the driving processes and the Xa's are 
idiosyncratic risk factors associated to the particular asset value. 

The basic Vasicek loan portfolio model assumes that the liabilities 
of the company are constant. However, this assumption is not always 
satisfied and liabilities may actually depend on the aggregate risk factor 
through the investment policy of the corporation (for an approach based 



on exponentially increasing but deterministic liabilities see Zhou (2001 1). 
Such a circumstance arises quite naturally, for example, in the case of 
unit linked insurance contracts. In order to introduce liability risk, we 
assume that the liabilities Bu evolve according to the random dynamics 
characterized by the stochastic differential equation 

dB it = ctiBudt + PiBit(V0ldYt + Vl - 6idZu), B i0 = B t , (3) 

where both the drift coefficient q 4 as well as the volatility coefficient j3i 
are exogenously given, 8i G [0, 1] is a coefficient measuring correlation 
between different liabilities, and Zit, ■ ■ ■ , Z„t are a family of independent 
driving Brownian motions independent of the aggregate risk factor Yt and 
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the asset-specific idiosyncratic risks Xu, . . . , X nt . Thus, the liabilities are 
assumed to fluctuate in a similar, yet not necessarily identical, fashion 
with the assets. 

As usually, we assume that default occurs whenever the assets do not 
meet the liabilities at a given date T. Since both the asset values as 
well as the liabilities follow two ordinary potentially correlated geometric 
Brownian motions, a standard application of Ito's lemma yields 

'A iT 



p [A iT < B iT ] = : 



< 1 



<3> 



IT 



where E? = of + 0} - 2ai0 V 1 OiQi measures the variance of the difference 
of the driving factors and 

S i = ln(g)-^-a i -i(<T?-^))r. 

For simplicity, we assume that the recovery rate from defaulted loans 
in the portfolio is zero. Under this assumption, the loss L; of the ith loan 
can be defined as the random variable 



Li = l(o,i) (AiT /Bir) 



1 if default occurs 
otherwise 

and, therefore, the loan portfolio percentage loss can be written as 

n jL - 1 

i — l 

Straightforward computation shows that the probability of default condi- 
tional on the aggregate factor Y now reads as 



p i (Y)=V[L i = l\Y]=\ 
where 



B, 



< 1 Y 



<1> 



A, 



Y 



C? = ff?(i-Pi) + ^(i-eo 

measures the variance of the difference of the idiosyncratic risk factors, 
and Aj = Oi^fpl — $i\[Wi denotes the volatility multiplier of the aggregate 
factor Y in the dynamics of the asset-liability ratio AiT /BiT- In contrast 
to the standard Vasicek loan portfolio model subject to deterministically 
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evolving liabilities, we now observe that the losses given default are inde- 
pendent random variables whenever the volatility multiplier A* is identi- 
cally zero for all the loans in the portfolio. As intuitively is clear, that 
case arises when the aggregate factor dynamics affects both assets as well 
as liabilities in a similar fashion. Otherwise, the losses are statistically 
dependent due to the joint dependence on the aggregate market factor. 
Moreover, applying the law of total probability shows that 

/oo p oo 

K[L i L j \Y T ]W[Y T €dy\ = / Pi (Vfy) Pj (VTy)$'(y)dy 
- oo J —oo 

implying that the covariance of the loss given default reads as 
cav[Li,L s ] = I" Pi (Vfy) Pj (Vfy)&(y)dy - $ ( -^=) $ 



Along the lines of our observations above, we find that if Ai = Aj = 
then cov[L,:, Lj] = 0. It is worth emphasizing that these covariances (and, 
therefore, default correlations) are typically very sensitive with respect to 
changes in the maturity T of the loans. 

In order to investigate the probability distribution of the percentage 
loss of a loan portfolio, let us now assume that the portfolio is formed by 
n identical contracts and denote the unconditional probability of default 
of an individual loan as p. In that case we observe that the probability of 
default conditional on the aggregate factor Y can be expressed as 

P(F) = $Q(£<£-- I (p)-AY)) ■ (4) 

Given this expression, denote now as L — lim n _>.oo L the limiting loan 
portfolio percentage loss of a infinitely large portfolio. We can now estab- 
lish the following: 

Proposition 2.1. The probability of k defaults in the loan portfolio per- 
centage loss reads as 



L — — 

n 



"» ' p k (VTy)(l-p(VTy)T- k $'(y)dy, (5) 



k 

where p(y) is given in Q . If A = then loan portfolio percentage loss 
converges almost certainly to the deterministic limit L — p. However, 
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i/A / then the limiting loan portfolio percentage loss is distributed 
according to the probability distribution 



F[L < x] = * (±- (C* -1 ^) - £<& _1 (p))) (6) 



with density 



/(a;) ~ |A| ' S'C*- 1 ^)) ' ' (7) 

Proof. The binomial formula ([5| is a direct implication of the law of total 
probability and the binomial nature of the loss given default (see, for ex- 



ample, chapter 9 in Lando (20041 and chapter 8 in |McNeil et al (20051). 
On the other hand, since the losses given default are conditionally inde- 
pendent, we observe that the conditions of the strong law of large numbers 
(SLLN) are satisfied and, therefore, that the percentage loss conditional 
on the aggregate factor converges to its expectation which, in the present 
case, reads as in Q when A / and as p when A = 0. Equation |6| then 
follows by computing the probability P[p(Y) < x]. The density can then 
be derived by ordinary differentiation. □ 



Proposition |2. 1| extends the results of the standard Vasicek loan port- 
folio model to the case where also liabilities are subject to random fluc- 
tuations. The main difference with the standard model is that now the 
volatility multiplier of the aggregate market factor in the dynamics of the 
asset-liability ratio can be zero even in the case where the factor affects 
both assets as well as liabilities. If this multiplier is zero, then the losses 
are IID random variables and the probability distribution can be directly 
analyzed in terms of constant binomial probabilities. In that case, the 
percentage loss converges almost everywhere to the known binomial prob- 
ability. However, if the multiplier is not zero, then the percentage loss 
converges towards a random variable with known distribution Q which 
resembles, but is not identical, with the limiting distribution in the case 
of constant liabilities. 
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Straightforward computations show that the probability density func- 
tion is bimodal when A 2 > £ 2 , monotone when A 2 = £ 2 , and unimodal 
with mode at 

^ = *(c^* -1 (P)) 
when A 2 < £ 2 . Consequently, along the original observations by Vasicek 
we find that depending on the precise parametrization of the model, the 
distribution may be either unimodal or bimodal and it can also be very 
skewed. We illustrate the loss density in Figure [l] for various correlations 
under the assumptions that 6 = 0.7, a = 0.2, /3 = 0.1, fj, = 0.055, a = 0.05, 
T = 1, Bo = 1, and A = 1.1. 




Figure 1: Loss densities 

In the present setting the ^-percentile L„ satisfying the identity f[L < 
L v ] = v is 

L^c&Q^-^ + IA^-V))). 

The percentile L v depends, among others, on the volatility fj of the liabil- 
ities. Unfortunately, it is not monotonic as a function of (3 and, therefore, 
the impact of liability risk on the percentiles is ambiguous. The 95% 
percentile Lo.95 is illustrated as a function of /3 in Figure [2] under the 
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assumptions 6 = 0.7, a = 0.2, fi = 0.055, a = 0.05, T = 1, Bo = 1, and 



A = 1.1. 




- p=0.7 

- p~0.83 

P =o.g 



03 - 



_J 1 1 L_ /J 

0.1 0.2 0.3 0.4 



Figure 2: The Impact of Liability Risk on the Percentile Lo.95 



3 The Effect of Random Systemic Risk 

Having considered the impact of liability risk on the limiting probability 
distribution of the loan portfolio percentage loss, we follow the original 



study by Zhou ( 19971 (see also Zhou (2001 1) and extend our basic model 
to the case where the assets backing up liabilities are subject to unex- 
pected random jumps modeled as a compound Poisson process. In con- 



trast with Zhou (19971, we assume that these unexpected jumps are only 
one-sided (downward jumps) and occur at the aggregate level. Therefore, 
the driving compound process is a common factor affecting all assets; an 
assumption permitting the analysis of the the impact of rare but poten- 
tially significant collapses (i.e. realization of systemic risk) in the asset 
values to the limiting default intensity in a large loan portfolio. 

In line with these arguments, we now assume that the asset values 
evolve according to the dynamics 



A, 



Ate 



(M l +A(l-E[ e -ei])-i CT f)t+ CTi (v^ y t + \/I^ x I t)-- 7 t 



(8) 
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where 

is a compound Poisson process independent of the continuous aggregate 
factor Y. In j9j, we assume Nt is a standard Poisson process with intensity 
A, {£fc}fc>i is a sequence of nonnegative iid random variables with known 
distribution, and £o = 0. In equation ([sjt A(l — E[e~^]) is a compensation 
term needed to guarantee that the asset value is expected to grow at 
the same rate as in the absence of jumps. If this compensation term is 
not taken into account then the proposed asset value model is almost 
surely lower and has a smaller expected value than the model considered 
in the previous section (due to the nonnegativity of the jumps and the 
monotonicity of the driving Poisson process). Especially, we observe that 
(Jsj) can be expressed as 

dA it = mAudt + (JiAudWit + A it / (e~ z - l)dN(dt, dz), (10) 

Jr 

where N(dt,dz) denotes the Poisson random measure associated to the 



underlying compensated Poisson process (cf. Chapter 2 in Kyprianou 



(2006)). 

It is worth pointing out that the stated specification results into an 
asset value which coincides in the mean but is more volatile than the 
model in the absence of systemic jumps. More precisely, it is clear that 
now that for all t it holds E[A it ] = Aie mi and 

vaxL4«] = A\e^ 1 f^+^nd-^if] _ ^ > A ^t l^t _ f 



In this way the considered process can be interpreted as a mean preserving 
spread of the continuous asset value dynamics considered in the previous 
section. 

Applying an analogous conditioning argument as in the previous sec- 
tion, we now find that the probability of default given the aggregate factors 
Y and J is 

¥[A tT < BiT | Y, J] = $ C_A= - ^=Y T 1 / '' 
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where Hi = S — A(l — E[e~^])T. As intuitively is clear, the positivity of 
the jump component Jt implies that the probability of default is in this 
setting higher that in the absence of unexpected downward jumps in the 
value of the assets. However, it is not beforehand clear how significant the 
effect of the Poisson component on the default probability is, and how this 
effect depends on both the intensity of the driving Poisson process and 
the precise nature of the jump size distribution. Moreover, in the present 
setting the loans are statistically dependent even when the volatility mul- 
tiplier of the aggregate market factor in the dynamics of the asset-liability 
ratio is zero (i.e. A; = for all i). The reason for this is naturally the 
presence of the systemic jump risk component affecting all assets. 

In order to be able to analyze the limiting probability distribution 
of the percentage portfolio loss, we now again assume that we have a 
portfolio of n approximately identical contracts. In that case we find that 
the conditional probability of default given the aggregate factors reads as 

^.fl-.(f (.-»>- + ^)), 



where 

/ S-A(l-E[e~ Sl ])T 
We can now establish the following result: 



p = <t> 



Proposition 3.1. If A = then the limiting loan portfolio percentage 
loss is distributed according to the probability distribution 

r{L<x] = e- XT Xlp .i ] (x) + J2 e ~ XT -T< L F[S k €du], (11) 

fc=i K - Jo 

where M T = £\/T (* _1 (a:) - $ _1 (p)) and 

' k 

P[S k G du] = 



.3 = 1 



= (g* ■ ■ ■ * g){u)du 



is the k-fold convolution of the density g(u) of the random jump-size. If, 
however, A / then 

P[L < x] = e~ XT ^(H(x,0)) 

+ Y. e ~ XT -^- $(H(x,u))F[S k £d U }, 
t-i fc - Jo 
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where 

/n this case, the density of the loan portfolio percentage loss reads as 
f (x) - _i_ Pfy^ (AT)* $>'(H(x,u)) . 

fix) -\A\J k] n*-H*)) [Sk€d] (13) 

— AT C 0)) 
|A| $'($-1(3;)) 

Proof. As in Proposition |2.1| the losses given default are conditionally 
independent given the aggregate factors and satisfy the conditions of 
the SLLN. The probability distributions |TT| and ( |T2"| follow directly 
by invoking the law of total probability in computing the probability 



P[p(y, J) < a;]. The density (131 can then be derived with ordinary dif- 



ferentiation. □ 



Proposition |3.1| states the limiting probability distribution and its den- 
sity for a sufficiently large loan portfolio percentage loss. Unfortunately, 
the distribution is in this case very complicated (being a mixture; for a 
comprehensive treatment of mixtures within credit risk management ap- 



plications, see Chapter 8 in McNeil et al (20051) and identifying the 



percentiles explicitly is extremely demanding, if possible at all. However, 
it is worth emphasizing that in contrast to the case subject to continuous 
factor dynamics, the distribution may now be multimodal. The reason 
for this observation is that now the density f(x) is a probability weighted 
sum of potentially bimodal densities. More precisely, since 

&(H(x,u)) 4*-i W a__i ?r ( f «-i W _ M -i CB _^) a 
$'($-!(a;)) 

is bimodal whenever A 2 > £ 2 , we notice that the limiting distribution may 
be multimodal depending on the jump size distribution. For example, 
when the jump size is a known constant, the limiting distribution may 
have more modes than just two. 

In order to investigate numerically the impact of jumps on the limiting 
distribution of the loan portfolio percentage loss, we now consider the spe- 
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cial case where the jump size is exponentially distributed with parameter 
7. It is well-known that in this case the series 



is Gamma-distributed according to the density 



P[S„ e du] = 



7e 



-du. 



(n-1)! 

In this case the density of the loan portfolio percentage loss reads as 
C r ^ .-at (AT) fc *'{H(x, u)) 7 e- 7tl (7^) fc - 



/(*) 



IAI 



fe! $'($-i( x )) (fe-1)! 



+ e 



-AT C ®'(H(X,0)) 



|A| $'($- 1 (a:))' 

We illustrate this density in the three different cases arising in the absence 
of jump risk. Figure [3] illustrates the case where the limiting distribution 
is unimodal under the assumptions that 9 = p — 0.7, a = 0.2, (3 = 0.1, 
iu = 0.055, a = 0.02, T = 1, A = 0.02, B = 1, and A) = 1.1. As is clear 




Figure 3: Loss densities 

from Figure [3j the presence of downward jump risk has a pronounced 
impact on the upper tail of the limiting density and, therefore, on the 
percentiles of the distribution. These percentiles are numerically illus- 
trated in the following table. As Table [l] clearly illustrates the difference 
between the percentiles is significant for sufficiently high percentiles. For 
example, in the absence of the systemic jump component the percentage 
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V 


0.90 


0.915 


0.93 


0.945 


0.96 


0.975 


7 — > oo 


57.1 


59.52 


62.23 


65.37 


69.12 


73.97 


7=1 


56.5 


59.32 


62.61 


66.6 


71.81 


80.01 


7 = 0.2 


54.65 


57.57 


61.01 


65.25 


70.98 


81.02 



Table 1: Percentiles in the Case of Figure [3] 



loss exceeds 73.97% with probability 2.5%. In the presence of the systemic 
jump component this percentile is radically changed and the percentage 
loss is expected to exceed 80.01% (81.02%) with the same probability. The 
expected shortfalls associated with the percentiles appearing on Table [I] 
are illustrated on Table [2] As Table [2] shows, the impact of the systemic 



V 


0.90 


0.915 


0.93 


0.945 


0.96 


0.975 


7 — > 00 


68.47 


70.26 


72.28 


74.61 


77.39 


80.97 


7 = 1 


72.7 


75.31 


78.39 


82.17 


87.09 


94.05 


7 = 0.2 


72.35 


75.22 


78.65 


82.9 


88.52 


96.43 



Table 2: Expected Shortfall ES„ 



jump component on the expected shortfalls is significant as well. Interest- 
ingly, the difference becomes higher as the confidence limit increases. The 
reason for this observation is the skewness of the density towards higher 
realizations in the presence of the systemic jump component. 

For the sake of comparison, the case where the limiting density is 
monotone in the absence of jump risk is illustrated in Figure [4] under the 
assumptions that 6 = 0.7, p w 0.83, a = 0.2, /3 = 0.1, fj, = 0.055, a = 0.05, 
T = 1, A = 0.02, B = 1, and A = 1.1. 

Again we notice from Table [3] that the impact of the discontinuous 
systemic risk component on the percentiles of the percentage loss distri- 
bution is significant. For example, in the absence of the systemic jump 
component the percentage loss exceeds 86.34% with probability 2.5%. In 
the presence of the systemic jump component these percentile is 91.69% 
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/to 




Figure 4: Loss densities 
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0.90 


0.915 


0.93 


0.945 


0.96 


0.975 


7 — > oo 


66.17 


69.42 


72.96 


76.85 


81.23 


86.34 


7=1 


65.94 


69.71 


73.91 


78.7 


84.38 


91.69 


7 = 0.2 


63.91 


67.89 


72.37 


77.59 


83.96 


92.84 



Table 3: Percentiles in the Case of Figure [4] 



(92.84%). The expected shortfalls associated with the percentiles appear- 
ing on Table [3] are now, in turn, illustrated on Table [3] 
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0.90 


0.915 


0.93 


0.945 


0.96 


0.975 


7 — > oo 


79.47 


81.54 


83.76 


86.18 


88.88 


91.98 


7=1 


82.26 


84.81 


87.61 


90.7 


94.18 


97.98 


7 = 0.2 


81.66 


84.45 


87.53 


90.97 


94.84 


98.91 



Table 4: Expected Shortfall ES V 



4 Conclusions 

We considered the impact of liability risk on the percentage loss distri- 
bution of a large uniform loan portfolio both in the presence and in the 
absence of discontinuous systemic risk. As our findings show, the impact of 
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liability risk is ambiguous and it may increase or decrease the percentiles 
depending on the precise parametrization of the considered model and, 
especially, on the strength of the dependence between the asset-liability- 
ratio and the the driving continuous aggregate factor. The discontinuous 
jump factor capturing the systemic risk has a more pronounced impact 
on the limiting percentage loss distribution since it affects all the asset- 
liability-ratios through the asset values. Our results seem to indicate that 
its impact becomes more significant at the tails of the distribution, which 
are found to be bimodal in the exponential case. According to our find- 
ings, the presence of systemic risk affects in a relatively significant way 
the expected shortfall associated to the upper tail probabilities even when 
the realization of the risk is assumed to be rare. 

There are several directions towards which our model could be gen- 
eralized. First, the considered loan portfolio is assumed to be large and 
uniform, thus overlooking the potentially significant effect of the gran- 
ularity of a loan portfolio. Second, assuming that there is no recovery 
once default has occurred is another simplifying assumption which could 
relaxed. Third, our analysis focuses solely on the distribution of the per- 
centage loss distribution and overlooks the pricing of bonds within the 
considered setting. All these interesting questions are left for future re- 
search. 
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